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ABSTRACT 


TJ 


The  effects  of  dynamic  perturbations  on  the  maximum  strength  of  initially 
imperfect,  axially  compressed,  circular  cylindrical  shells  are  studied 
using  a  modified  Hamilton's  principle  based  on  a  Reissner-type  functional, 
von  Karman-Donnell  shall  kinematics,  and  a  deformation  theory  of  plas¬ 
ticity.  The  dynamic  effects  have  been  included  by  making  use  of  the 
analogy  of  time-dependent  imperfection  growth  to  represent  dynamic 
disturbance  of  the  lateral  motion  of  the  shell  wall.  The  results  of 
the  present  analysis  reflect  families  of  load  versus  end-shortening 
curves  for  long  circular  cylinders.  For  each  material,  initial  imper¬ 
fection  parameter,  and  radius-to-wall  thickness  ratio,  there  exist 
frequency-shortening  rate  ratios  which  provide  unique  load-end  short¬ 
ening  paths  in  both  the  elastic  and  inelastic  ranges.  Significant 
maximum  strength  reductions  are  obtained  for  aluminum  and  stainless 
steel  relative  to  the  predictions  of  static  elastic  and  inelastic 
theoretical  analyses  for  shell  radius-to-thickness  ratios  represen¬ 
tative  of  closely  stiffened  and  sandwich  construction. 
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FOREWORD 


The  work  reported  herein  constitutes  a  portion  of  the  effort  being  un¬ 
dertaken  at  Stanford  University  for  the  Eustis  Directorate,  U.S.  Army 
Air  Mobility  Research  and  Development  Laboratory  under  Contract 
DAAJ02-70-C-0075  (Task  1F061102A33F  02’) .  The  present  program  is 
one  of  a  series  aimed  at  establishing  accurate  theoretical  predic¬ 
tion  capability  for  the  static  and  dynamic  behavior  of  aircraft 
structural  components  using  both  conventional  and  unconventional 
materials.  Predecessor  contracts  supported  investigations  which 
led,  in  part,  to  the  results  presented  in  References  10,  13,  14, 

15,  20,  38,  and  31. 
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INTRODUCTION 


The  theoretical  and  experimental  developments  pertaining  to  thin  shells 
in  recent  years  have  essentially  been  channeled  into  three  distinct 
areas:  buckling  and  postbuckling,  free  vibration,  and  impulsive  or 
forced  excitation.  In  each  case,  attempts  have  been  made  to  resolve 
the  discrepancies  between  classical  theory  predictions  and  experimental 
results.  In  general,  developments  in  each  area  have  taken  place  inde¬ 
pendently  of  the  others. 

Much  of  the  work  directed  toward  initial  buckling  and  postbuckling 

12  3 

has  been  surveyed  in  three  comprehensive  papers  by  Hoff  ’  and  Stein. 

The  strongest  emphasis  by  investigators  has  been  directed  toward  asym¬ 
metric  buckling  of  thin-walled,  circular  cylindrical  shells  in  axial 
compression.  The  consensus  of  both  authors  has  attributed  the  dis¬ 
crepancies  between  theory  and  experiment  for  compressed  shells  to  the 

combined  effects  of  initial  imperfections,  prebuckling  deformations,  and 

4 

boundary  conditions.  The  theoretical  work  of  Donnell  and  Wan  and 

Koiter,^  for  example,  coupled  with  the  experimental  findings  of  Horton 

6  7  8 

and  Durham,  Babcock  and  Sechler,  and  Almroth,  Holmes  and  Brush0  has 

amply  demonstrated  the  significance  of  the  presence  of  initial  imper¬ 
fections  in  the  buckling  of  axially  compressed  shell  structures.  Lesser 
significance  is  attributed  to  prebuckling  deformations  and  boundary  con¬ 
ditions  from  either  quantitative  or  practical  considerations.  However, 
not  to  be  overlooked  is  the  effect  of  inelastic  deformations;  this 
phenomenon  was  apparently  appreciated  by  Donnell^  in  195^  but  not  given 
serious  attention  until  the  past  decade,  commencing  with  the  theoretical 
work  of  Mayers  and  Rehfield^  and  the  experimental  program  of  Edwards. ^ 

With  the  use  of  variational  principles,  nonlinear  stress-strain  relations, 

and  von  Karman's  large-deflection  strain-displacement  relations,  Mayers 

12  13 

et  al.  ’  were  able  to  show  for  axially  compressed  flat  rectangular 
plates  that  unique  material  dependent  and  width-to-thickness  dependent 
load-end  shortening  curves  existed  when  compression  beyond  the  classical 


r 


buckling  load  and  into  the  inelastic  range  occurred.  More  recently, 

14  15 

Mayers  and  Wesenberg  *  have  been  able  to  obtain  material-dependent 
and  radius-to-thickness  ratio-dependent  load-shortening  curves  for 
axially  compressed  cylinders  with  the  use  of  a  modified-Reissner  var¬ 
iational  principle,  von  Karman- Donne 11  strain-displacement  relations, 
and  nonlinear  constitutive  relationships.  Unlike  their  plate  counter¬ 
parts,  the  families  of  initial-imperfection-dependent  load-shortening 
curves  never,  in  reality,  reach  classical  buckling  load  predictions. 

In  the  presence  of  imperfection-sensitive  loading  paths  into  the  in¬ 
elastic  range,  significant  reductions  in  maximum  strength  for  both 
stiffened  and  unstiffened  shells  occur.  Such  reductions  are  shown 
to  be  obtained  relative  to  any  theory  making  use  of  a  linearized 
stress-strain  law.  Significant  reductions  result  even  though  the 
average  compressive  stress  at  buckling  is  well  below  the  0.2$  offset 
yield  stress  of  the  material  and  as  low  as  the  proportional  limit 
stress  for  some  materials. 

With  regard  to  vibration  of  thin  shells,  much  of  the  early  development 
done  within  the  context  of  linearized  theory  is  contained  in  the  work 
of  Rayleigh,  ^  Love,*^  and  Flu'gge.^  In  1955,  Reissner^  made  the  first 
significant  contribution  to  the  solution  of  the  nonlinear  problem  by 
showing  the  axial  and  circumferential  inertia  to  be  negligible  for 
predominantly  radial  motion.  It  is  noteworthy  to  mention  that  his 
results  simplified  to  those  derivable  from  von  Karman -Donne 11  theory. 
Reissner,  however,  selected  the  asymmetric  chessboard  deflection  pat¬ 
tern  of  linear  theory  for  his  shell  nonlinear  vibration  studies.  This 

choice  has  been  invalidated  by  subsequent  work  on  periodic,  nonlinear 

„  20 

free  vibrations  of  thin  shells.  Finally,  Mayers  and  Wrenn  have  shown 

that  nonlinear,  free  vibrations  of  shells  are  nonperiodic  in  nature 

based  upon  minimum  energy  considerations. 

The  direction  of  this  study  has  been  to  combine  the  effects  of  static 
maximum  strength  and  vibration  in  the  more  encompassing  area  of  dynamic 
stability.  A  number  of  subclasses  of  the  field  of  dynamic  stability 


2 


are  mentioned  here.  Each  combines  the  effects  of  external  loading  and 

kinematic  displacement  of  the  structure  to  obtain  a  particular  effect. 

21 

Varying  combinations,  as  noted  by  Hoff,  have  been  classified  as  para¬ 
metric  resonance,  impulsive  loading,  circulatory  loading,  aeroelasticity, 
and— the  one  most  pertinent  to  this  study— buckling  in  the  test  machine. 
The  phenomena  of  resonance  and  aeroelasticity  have  been  well  studied 
through  the  years.  They  are  related  to  the  excitation  of  one  of  the  nat¬ 
ural  frequencies  of  the  s tructure  in  ques  tion  until  such  time  as  the  ampli¬ 
tude  range  surpasses  some  admissible  bound.  Impulsive  loadings,  on  the 
other  hand,  have  traditionally  been  associated  with  the  massive  applica¬ 
tion  of  load  for  short  durations  beyond  the  maximum  permissible.  Cir¬ 
culatory  loading  is  not  considered  here  due  to  the  presence  of  special 
highly  directional  loading  requirements. 

Dynamic  buckling  criteria  in  an  elastic  testing  machine  have  been  studied 

22 

in  several  specialized  cases.  Hoff  has  shown  the  existerce  of  constant 
total  energy  contour  lines  about  three  equilibrium  points  for  a  com¬ 
pressed  imperfect  column.  Similarly,  by  approaching  the  problem  through 

consideration  of  the  effects  of  loading  velocity  in  a  testing  machine, 
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Hoff  was  able  to  detect  stress  reversal  before  maximum  load  in  the 
case  of  columns. 

In  the  present  work,  the  maximum  strength  behavior  of  axially  compressed 
circular  cylindrical  shells  in  the  presence  of  dynamic  perturbations  has 
been  studied  using  a  modified  Hamilton's  principle  incorporating  time- 
dependent  imperfection  growth  as  a  suitable  analog.  Also  incorporated 
as  part  of  the  variational  principle  is  a  modified  form  of  Reissner's 
functional  based  upon  the  von  Karman-Donnell  strain  displacement  rela¬ 
tions.  The  allowance  for  the  presence  of  the  modified-Reissner  functional 
facilitates  the  incorporation  of  inelastic  effects  by  selection  of  the 
state  of  stress  independently  of  the  displacements  (or  state  of  strain). 
This  procedure  has  been  well  established^  and  used  with  confidence  for 
shells  inReferences  l4andl5.  On  this  basis,  inelastic  maximum  strength, 
radius-to-thickness  dependent  load-shortening  curves  have  been  obtained. 
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The  effects  of  the  presence  of  dynamic  perturbations  are  then  studied 
not  only  in  terms  of  the  mechanism  of  buckling  but  also  in  terms  of  both 
shell  radius-to-thickness  ratio  and  material  property  dependency. 
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GENERAL  THEORY 


STATEMENT  OF  PROBLEM  AND  BASIC  ASSUMPTIONS 


The  problem  considered  is  that  of  the  effect  of  dynamic  perturbations 
in  the  maximum  strength  analysis  of  initially  imperfect,  axially  com¬ 
pressed  circular  cylindrical  shells.  The  end-shortening  rather  than 
the  end-loading  of  the  shell  is  prescribed;  this  corresponds  to  loading 
of  the  shell  in  a  rigid  testing  machine.  The  approach  to  the  solution 
has  been  to  consider  the  physical  interaction  of  free  vibration  in  the 
deflected  shape  of  the  cylinder  wall  as  a  means  of  initiating  growth 
in  the  initial  imperfection  shape.  The  differential  equations  of  mo¬ 
tion  are  reduced  to  their  static  equivalents  and  an  additional,  approx¬ 
imate  sinusoidally  vibrating  deflection  function  is  included  to  account 
for  the  effect  of  inertia. 

The  solution  for  the  differential  equations  of  motion  for  axially  com¬ 
pressed,  vibrating  shells  is  obtained  from  a  modified  Hamilton's  vari¬ 
ational  principle.  Reissner's  functional  is  substituted  for  the  strain 
energy  in  order  to  provide  arbitrariness  in  the  selection  of  strain  and 
stress  states  and  to  include  the  effects  of  nonlinear-elastic  material 
behavior.  In  the  statically  equivalent  system,  no  unloading  in  the 
nonlinear-elastic  range  is  assumed  to  occur.  In  the  assumed  dynamic 
deflection  function,  only  the  first  quarter-sine  wave  is  considered. 
Since  maximum  strength  is  the  criterion  for  comparison,  the  first 
maximum  amplitude  is  all  that  need  be  included.  As  a  result,  the  in¬ 
elastic  and  nonlinear-elastic  ranges  can  be  considered  as  interchange¬ 
able.  From  the  reduction  of  the  differential  equations  of  motion  to 
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their  static  equivalents,  the  results  of  Mayers  and  Wesenberg  are 
reproduced. 

The  basic  two-element  description  of  the  shell  cross  section  (Figure  l) 
previously  used  in  plate  and  shell  analyses  (see  References  12,  13,  14, 
15)  is  adopted  for  the  dynamic  maximum  strength  analysis.  This  descrip¬ 
tion  in  its  simplest  form  is  composed  of  two  thin  face  sheets  separated 
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Figure  1 


Circular  Cylindrical  Shell  With 
Two-Element  Cross  Section. 
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by  a  non-stress-carrying  core.  The  stress  distribution  is  considered  to 
be  constant  across  the  thickness  o f  any  face  sheet,  and  the  complexities 
of  integrating  a  nonlinear  stress  distribution  through  the  shell  thick¬ 
ness  are  conveniently  avoided.  The  core  is  considered  to  be  infinitely 
rigid  in  shear  and  normal  compression;  the  shell  behavior  is  based 
on  von  Karman-Donnell  large  displacement  theory.  In  general,  ma¬ 
terial  compressibility  effects  are  considered  to  be  negligible  and 
Poisson's  ratio  is  taken  as  v  =  0.5  for  computations.  The  constitu¬ 
tive  equations  used  to  describe  the  nonlinear  material  properties  make 
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use  of  the  Ramberg-Osgood  power  law  stress-strain  curve  approximation 
for  aircraft  structural  metals. 

The  von  Karman-Donnell  strain-displacement  equations  have  been  adopted 
for  the  two-element,  initially  imperfect  shell.  For  radial  deflection- 
thickness  ratio  magnitudes  of  the  order  of  one  in  the  prebuckling  or 
initial  buckling  state,  and  in  view  of  the  conclusions  of  Mayers  and 
Rehf ield, '*'0  the  von  Karman-Donnell  strain-displacement  relations  are 
justifiable. 


BASIC  EQUATIONS 


The  von  Karman-Donnell  strain-displacement  relations  for  the  isotropic 
two-element  cylinder,  modified  to  include  the  effects  of  initial  radial 
imperfections,  as  shown  in  Reference  lh,  are 
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The  curvature-displacement  relations  are  given  by 
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The  total  strains  for  the  top  and  bottom  faces  of  the  two-element  iso¬ 
tropic  cylinder  become,  respectively, 
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REISSNER  FUNCTIONAL 


(3) 


The  Reissner  variational  principle  requires  the  vanishing  of  the  first 
variation  of  both  the  free  states  of  stress  and  strain  for  equilibrium 
and  stress-displacement  compatibility.  The  requirements  for  admissible 
variations  are  the  enforcement  of  the  boundary  conditions  on  the  pre¬ 
scribed  surface  displacements  and  their  continuity  throughout  the  body. 
For  plates  and  shells,  the  Reissner  functional  is 


u"  =  fff  Jcr  e  +  a  e  +  t  y  -F 
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where  F'  is  the  stress-energy  density.  F'  is  the  energy  function 
which  for  incremental  changes  in  stress  relates  the  strains  at  every 
point  by  the  expressions 


dF' 

Lx  cV 


dF^ 

cVr  ’ 

y 


xy 


(3) 


8 


L 


f 

KINETIC  ENERGY 

The  two-dimensional  kinetic  energy  expression  is  the  sum  of  the  incre¬ 
mental  longitudinal,  circumferential,  and  lateral  motions  of  the  cylin¬ 
der.  For  the  analysis  of  thin  shells,  rotational  inertia  effects  are 
neglected.  The  kinetic  energy  then  becomes  the  sum  of  the  kinetic 
energies  associated  with  the  axial,  circumferential,  and  radial 
velocities  of  the  cylinder,  or 

T"  =  J£f  2  p  u2  +  v2  +  w2  dxdydz  (6) 

V 

Integration  over  the  depth  of  the  two-element  cylinder  of  face  thick¬ 
nesses  tf  gives 

L  2rrR 

t"  =  |  p(2tf)  J  J  [u2  +  V2  +  V2  dxdy  (7) 

0  0 

for  constant  p 
STRESS-ENERGY  DENSITY 

For  linear-elastic  problems,  F*  becomes  the  complementary  energy 
density  and  Reissner's  functional  can  be  reduced  to  the  strain  energy 
with  the  enforcement  of  Hooke's  law.  As  developed  in  Reference  12 
for  nonlinear  elastic  materials,  the  stress  energy  density  F/  is 
•»iven  by 

°eff 

F'  =  /  eeffdtTeff  (8) 

0 

where  from  the  secant-modulus  deformation  theory  of  plasticity  the 
"effective"  stress  is  written  as 

cr2  £  =  a2  +  cr2  -  2w  cr  +  2(1  +  v)  (9) 

eff  x  y  v  x  y  v  '  xy 

for  an  isotropic  material  and  v  is  understood  to  equal  \  . 
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In  turn,  "effective"  stress  is  related  to  "effective"  strain  by  the 
secant-modulus  of  the  uniaxial  stress-strain  curve  of  the  material  as 


°ef£  ESeeff 

In  the  present  analysis  developments,  the  uniaxial  stress-strain  curves 
in  tension  and  compression  are  taken  as  identical. 


Upon  substitution  for  the  stress  energy  density,  Reissner's  functional 
may  then  be  written  as 
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For  the  purposes  of  this  analysis  the  nonlinear  relationship  between 
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stress  and  strain  is  that  developed  by  Ramberg  and  Osgood.  This 
three-parameter  representation  of  the  uniaxial  stress-strain  curve  is 
written  as  a  power  function  in  the  form 


where  E  ,  K  ,  and  N  are  the  material  constants  for  a  given  material. 
For  the  present  analysis,  this  form  is  used  without  regard  to  the  ef¬ 
fects  of  unloading.  Previous  work  involving  inelastic  effects  in  the 
bending,  buckling,  and  postbuckling  of  plates  and  shells  indicates  that 
only  localized  unloading  of  small  magnitude  takes  place  at  buckling. 


The  functional  u"  integrates  to 
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Consistent  with  the  two-element  description  developed  in  Appendix  I, 
integration  of  Reissner's  functional  over  the  depth,  h  ,  yields 
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in  terms  of  the  midsurface  direct  and  bending  stress,  the  strain  com¬ 
ponents,  and  the  “effective"  stress 


=  2t. 
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where  the  "t"  and  "b"  subscripts  refer  to  the  top  and  bottom  faces 
of  the  two-element  cross  section,  respectively.  Finally,  with  the 
addition  of  the  von  Karman-Donnell  strain-displacement  relations,  the 
Reissner  functional  becomes 
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The  kinetic  energy  is  combined  with  the  Reissner  functional  in  the  form 
of  a  modified  Lagrangian.  When  integrated  over  a  given  time  period, 
this  yields  the  Hamiltonian  in  the  form 
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VARIATIONAL  PRINCIPLE 


The  application  of  the  mod  if ied-Hamil ton's  principle  requires  the 
vanishing  of  the  first  variation  simultaneously  with  respect  to  the 
admissible  degrees  of  freedom  that  characterize  the  states  of  stress 
and  strain;  this  condition  is  indicated  by  the  symbolic  form 
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(U"  -  T")  dt  =  0 


(17) 


The  prescribed,  constant-rate  end  shortening  e(t)  is  included  in 
u(t)  ;  thus,  a  potential  of  applied  loads  term  does  not  appear  in 
the  variational  principle. 


The  application  of  this  principle  allows  the  independent  variation  of 
the  states  of  stress  and  strain.  The  enforcement  of  the  above  condi¬ 
tion  results  in  the  Euler  equations  of  motion,  the  stress-displacement 
relations,  and  both  the  initial  and  boundary  conditions.  The  system 
of  Euler  equations  for  the  circular  cylindrical  shell  which  behaves 
according  to  the  von  Karman-Donnell  strain-displacement  relations  is 
developed  in  Appendix  I  for  the  two-element  model  representation  of 
the  actual  cylinder. 
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METHOD  OF  SOLUTION 


Investigators  in  the  past  have  approached  the  maximum  strength  and 

nonlinear  vibration  problems  as  two  distinct  phenomena.  In  References 

1-3,  the  major  contributors  to  the  linear-elastic  buckling  and  post- 

buckling  of  shells  to  1967  are  listed.  More  recently,  Mayers  and 

14  15 

Wesenberg  ’  have  extended  theory  and  analysis  to  include  inelastic 
effects.  For  nonlinear,  free  vibrations  of  shells,  the  literature  is 

relatively  sparse.  Noteworthy  contributions,  all  since  1955>  are  those 

19  2k  20 

of  Reissner,  Evenson  and  Fulton,  and  Mayers  and  Wrenn. 

The  present  analysis  has  taken  the  perspective  of  treating  vibration 

as  a  time-dependent  perturbation  of  initial  imperfection  in  a  shell 

which  decreases  the  maximum  strength.  Thus,  the  analysis  can  proceed 

14 

much  like  that  of  Mayers  and  Wesenberg  in  its  development  of  a  so¬ 
lution.  A  prescribed  constant-rate  displacement  of  the  shell  edges 
permits  an  assumed  time-dependent  imperfection  growth  (perturbation) 
amplitude  to  be  related  to  the  end  shortening. 

DISPLACEMENT  FORMULATIONS 


The  primary  consideration  in  the  choice  of  a  displacement  function  is 
the  prebuckling  deformation  pattern.  The  approach  followed  is  to  take 

the  free  stresses  and  radial  deflection  pattern  in  the  same  form  as 

25 

those  formulated  by  Kempner.  Thus,  the  oresent  analysis  can  be  re¬ 
duced  to  that  of  Kempner' s  in  the  perfect  cylinder,  linear-elastic, 
static  case.  Then,  rather  than  specifying  the  stresses  and  strains 
as  given  functions  of  time,  the  unknown  coefficients  are  allowed  to 
vary  freely  and  the  differential  equations  of  motion  in  terms  of  the 
initial  and  buckling  patterns  are  obtained. 

25 

The  radial  displacement  buckle  pattern  as  specified  by  Kempner  is 


w(t)  =  h 
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where  the  wavelengths  are  X  =  L/m  and  X  =  ?rR/n  ,  respectively. 

x  y 

For  prescribed  end  shortening,  the  choice  of  the  corresponding  midsur¬ 
face  displacements  is 


u  =  -  ex 

v  =  0  (19) 

where  e  =  et  is  the  unit  end  shortening  and  e  is  the  constant  unit 
end-shortening  rate. 


The  assumption  of  midsurface-displacement  functions  reflecting  no  free 

parameters  has  been  justified  in  References  10,  13,  14,  and  15.  This 

simplification  takes  cognizance  of  the  fact  that  either  nonexistent  or 

only  very  weak  coupling  can  be  assumed  to  occur  between  the  midsurface 

stresses  and  strains  in  the  Reis^ner  formulation  of  the  shell  bending 

and  buckling  problem  with  only  lateral  equilibrium  being  satisfied  by 

application  of  the  variational  principle.  An  analogous  simplification 
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in  the  shell  vibration  problem  is  attributed  to  Reissner;  as  long  as 
the  inertia  effects  are  primarily  the  results  of  transverse  vibration, 
longitudinal  inertia  effects  can  be  neglected.  As  shown  in  Appendix  I, 
the  simplification  in  the  u-  and  v-displacement  choices  permits 
satisfaction  of  the  middle-surface  equilibrium  equations  independently 
of  the  magnitudes  of  the  free  stress  coefficients. 


STRESS  FORMULATIONS 


The  membrane  stresses  are 
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As  shown  in  Appendix  I,  this  set  of  stresses  identically  satisfies  the 
middle-surface  equilibrium  equations  for  the  x-  and  y-directions . 


The  bending  stresses  are  chosen  in  the  form  of  the  curvatures,  but 
with  free  coefficients  to  accommodate  the  general  case  of  a  nonlinear 
relationship  between  stress  and  strain.  Thus, 
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INITIAL  IMPERFECTION  PATTERN 

26 

In  the  experimental  work  of  Tennyson  and  Welles,  the  first  approxima¬ 
tion  to  the  radial  deflection  pattern  at  the  onset  of  buckling  was  found 
to  be  composed  of  the  asymmetric  chessboard  pattern  plus  an  axisymmetric 

contribution  corresponding  to  the  second  harmonic.  Similarly,  and  with 

20 

more  significance  to  the  present  work,  Mayers  and  Wrenn,  in  their 
analysis  of  the  nonlinear  free  vibrations  of  thin  shells,  determined 
that  for  radial  deflections  of  the  order  of  the  shell  thickness,  the 
minimum  energy  solution  corresponds  to  the  chessboard  pattern  modified 
by  a  symmetric  contribution  corresponding  to  the  second  harmonic. 

The  choice  of  the  imperfection  pattern  herein  includes  two  components: 

26 

the  initial  buckling  shape  ascribed  to  Tennyson  and  Welles  and  used 

14 

by  Mayers  and  Wesenberg  in  their  maximum  strength  analysis,  and  an 
imperfect  ion -growth  component  with  the  sane  modal  shape  but  with  time- 
varying  amplitude.  It  is  this  prescribed  growth  component  that  provides 
the  analogy  to  the  presence  of  inertia  in  the  equations  of  motion.  The 
form  of  the  initial  imperfection  is 


where,  as  a  first  approximation  to  dynamic  perturbation  on  static 
maximum  strength,  simple  harmonic  motion  is  assumed  to  represent  the 

imperfection  growth  pattern.  Though  nonlinear  vibration  of  circular 
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shells  has  been  shown  to  be  nonperiodic,  the  primary  interest  here 
is  in  carrying  the  analysis  through  the  first  quarter-period  of  an 
imperfection-growth  oscillation.  At  this  time,  the  perturbation  has 
reached  its  first  maximum  amplitude. 
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ELASTIC  LOAD-SHORTENING  F^LATIONSHIP 

With  the  displacement,  stress,  and  imperfection  functions  given  by 
Equations  (l8),  (19),  (20),  (21),  and  (22),  all  of  the  quantities 
appearing  in  the  functional,  Equation  (l6),  are  defined  and  the 
variational  principle  indicated  by  Equation  (17)  can  be  applied. 

The  integrated  functional  is  given  in  Appendix  II  for  the  general 
case.  With  time-dependence  neglected  (static  problem)  and  K  =  0 
(linear-elastic  material),  the  Reissner  functional  reduces  to  that 
developed  originally  in  Reference  14  and  reproduced  in  Appendix  III. 

For  the  present  static  problem,  with  K  =  0,  the  vanishing  of  the 
first  variation  of  Equation  (l6)  with  respect  to  |  ,  , 

bij,  dj^  ,  and  cr/E  leads,  after  simplification,  to  just  four  equations 
in  the  four  unknowns,  |  ,  and  cr/E  .  Their  solution 

for  given  imperfection  growths  nnd  end-shortenings  permits  the  construc¬ 
tion  of  the  linear-elastic  material  load-shortening  curves  shown  in 
Figure  2. 

INELASTIC  LOAD-SHORTENING  RELATIONSHIP 


The  equivalent  inelastic  static  problem  is  obtained  for  nonvanishing 
K  (nonlinear  elastic  material)  in  Equation  (l6).  The  inelastic  load' 
shortening  curves  for  the  circular  cylindrical  shell  are  derived  from 
the  extrema  of  the  following  functional: 


U" 

(-) 

1  K  11 

-  If 

r  N+l  v  N+ln 

(^eff\  +/-S«) 

VE 

^  elastic 

2  (N  +  1 )  { J0 

L  V  E  /  t  \  E  /  b 

•*(:)•(=) 

The  vanishing  of  the  first  variation  with  respect  to  each  of  the  17 
free  parameters  ,  A„  ,  a^  ,  b^  ,  d^  ,  and  a/E  is  enforced 

to  again  obtain  a  system  of  imperfection-dependent  load-shortening 
relationships.  The  presence  of  higher-order  terms  associated  with 

18 


L 


<rR/Et 


nonzero  values  of  K  precludes  any  rational  possibility  of  obtaining 
explicit  relationships  coupling  any  of  the  17  free  parameters.  The 
solution  for  the  extrema  has  been  adopted  for  numerical  solution  using 
a  Newton-Raphson  technique  with  double  integration  over  the  surface  of 
any  one  buckle  of  the  shell.  The  resulting  maximum- strength  radius-to- 
thickness  dependent  load-shortening  curves  have  been  obtained  for  2024-T3 
aluminum  and  stainless  steel  (^-  hard),  the  stress -strain  curves  for  which 
are  shown  in  Figures  3  and  4,  respectively.  The  load-shortening  curves 
are  presented  in  Figures  5,  6,  7,  8. 

THE  DYNAMIC  ANALOGY  -  IMPERFECTION  GROWTH 

Imperfection  growth  is  achieved  by  the  superposition  of  time-dependent 
imperfection  components  on  the  initial  static  imperfection  pattern 
components.  As  indicated  by  Equation  (22),  simple  harmonic  motion 
is  selected  as  describing  the  behavior  of  the  dynamic  disturbance 
occurring  in  the  initially  imperfect,  uniformly  shortened  shell. 

The  linear  combination  of  the  static  and  dynamic  imperfections  be¬ 
comes  the  total  imperfection  at  any  instant  in  time.  In  view  of 
Equations  (19)  and  (22),  the  disturbance  or  imperfection  growth 
and  the  unit  end  shortening  e  are  coupled  in  time  to  create  a 
phenomenon  analogous  to  end-loading-excited  lateral  vibration. 

With  this  analog,  the  effect  of  small  dynamic  disturbances  on  the 
static  maximum  strength  of  initially  imperfect  shells  can  be  studied. 

The  manner  of  utilization  of  the  analog  ir.  the  present  investigation 
is  illustrated  in  Appendix  IV. 
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in  Present  Inelastic  Analysis. 
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Figure  6.  Family  of  Static  and  Dynamic  Imperfection-Sensiti 
Load-Shortening  Curves  for  a  Two-Element,  202U-T3 
Aluminum  Cylinder  With  R/t  =  100. 
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DISCUSSION  AND  RESULTS 


The  maximum  strength  of  dynamically  perturbed,  initially  imperfect, 
axially  compressed,  circular  cylindrical  shells  has  been  studied 
using  an  analog  to  represent  dynamic  imperfection  growth.  A  modi¬ 
fied  Hamiltonian  based  on  the  von  Karman- Donne 11  strain-displacement 
relations  and  a  nonlinear  constitutive  law  is  incorporated  to  obtain 
Euler  equations  of  motion  of  the  nonlinear,  though  conservative,  sys¬ 
tem.  The  effects  of  inertia  in  the  equations  of  motion  are  neglected, 
and  a  dynamic  imperfection  component  is  added  to  the  static  initial 
imperfection  to  obtain  a  simple  analog  of  the  dynamically  perturbed 
problem.  Maximum  strength,  radius-to-thickness  dependent,  load¬ 
shortening  curves  are  obtained  for  both  the  static  problem  and  the 
analog  of  the  dynamic  problem.  Two  different  radius-to-thickness 
ratios,  100  and  200,  are  used  to  demonstrate  the  effects  on  maximum 
strength  in  the  elasto-plastic  region.  Similarly,  two  materials, 
stainless  steel  hard)  and  2024-T3  aluminum  alloy,  are  used  as 
examples  of  practical  materials  with  vastly  different  nonlinear 
constitutive  relations.  Several  different  imperfection  ranges  are 
selected  to  correlate  qualitatively  the  results  of  experimental  work 
undertaken  during  the  past  several  decades. 

ELASTIC  LOAD-SHORTENING  RELATION 


A  family  of  imperfection-sensitive,  static,  linear-elastic,  load¬ 
shortening  curves  for  an  isotropic,  circular  cylindrical  shell  have 
been  produced  in  Figure  6  utilizing  von  Karman-Donnell  large- 
displacement  theory  in  a  modified-Reissner  variational  process. 

A  Rayleigh-Ritz  procedure  based  on  a  four-term  trigonometric  series 
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of  the  same  form  as  that  of  Kempner,  has  been  used  in  conjunction 
with  the  linearized  moment-curvature  relations  and  the  middle-surface 
equilibrium  equations  to  obtain  the  linear-elastic  solution  of  the 
general  problem.  The  results  obtained  represent  the  solution  to  the 
initially  imperfect  problem  presented  in  Appendix  II.  The  values  of 
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the  waveform  parameters,  4  =  1.0  and  T)  =  0.826,  are  selected  as 

stationary  values  for  the  bending  and  buckling  range  based  on  solu- 
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tion  results  of  (l)  the  linearized  problem  by  Madsen  and  Hoff  and 
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(2)  the  nonlinear  problem  by  Meller  and  Mayers,  who  obtained  these 
values  in  the  limit  as  the  bifurcation  point  of  linear  theory  is 
approached  for  a  perfect  shell.  These  values  are  maintained  during 
the  imperfect  shell  analysis  in  view  of  their  collective  minimum- 
energy  property  for  the  perfect  shell. 

The  reduction  in  maximum  strength  with  increasing  initial  imperfection 
is  clearly  shown  in  Figure  2.  As  small  an  initial  deviation  of  1$  of 
the  wall  thickness  from  the  perfect  shape  can  produce  a  15$  reduction 
in  Initial  buckling  load.  For  an  Initial  imperfection  parameter  value 
of  about  ^11q  =  0.05,  the  load-shortening  curve  changes  from  a 
multivalued  relationship  at  constant  eR/t,  in  which  catastrophic 
snapthrough  is  the  characteristic  failure  mode,  to  a  single-valued 
function  at  constant  eR/t  with  rapid,  but  continuous,  reduction 
in  load  versus  end  shortening,  and  about  a  3C$  reduction  in  buckling 

strength.  The  reduction  in  buckling  strength  continues  until,  as 

28 

shown  by  Mayers  and  Meller,  the  concept  of  initial  maxima  vanishes 

completely  for  ^11  ~  0.17.  It  must  be  noted  that  the  use  of  the  four 
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term  deflection  pattern  of  Kempner  produces  accurate  bending,  and 
initial  buckling  results;  but  once  the  initial  maxima  have  been  reached 
the  limitations  of  von  Karman-Donnell  theory  are  surpassed  and  the 
solutlo  fails.  It  was  shown  first  by  Yoshimura  that  the  very  large- 
deflection  postbuckled  shape  of  the  circular  cylindrical  shell  is  de¬ 
velopable  (no  stretching  of  the  middle  surface)  into  a  polyhedral  sur¬ 
face  consisting  of  Identical  plane  triangles.  In  a  study  of  the 
validity  of  von  Karman-Donnell  shell  theory,  Mayers  and  Rehfield^ 
showed  the  theory  to  be  completely  inadequate  for  representing  very 
large -deflection  postbuckling  behavior  except  in  the  limiting  case 
of  linear  elasticity  and  t/R  -•  0  . 
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In  the  survey  papers  of  Hoff  ’  and  Stein,  thorough  review  is  made  of 
the  possible  reasons  for  the  significant  discrepancies  present  between 
theory  and  experiment  in  shell  buckling  under  axial  compression.  The 
scatter  band  of  test  data  reported  through  1961  (see  Reference  3)  is 
shown  in  Figure  9-  To  obtain  physically  representative  results  in 
the  present  analysis,  a  range  of  imperfection  values  based  on  linear- 
elastic-analysis  calculations  has  been  selected  to  yield  typical 
upper,  middle,  and  lower  bandwidth  values  of  shell  maximum  strengths. 

A  2024-T3  aluminum  alloy  cylinder  with  R/t  =  200  is  used  as  a 
representative  elastic  shell,  since  its  buckling  stresses  based  on 
inelastic  analysis  are  within  99$  of  the  linear-elastic-analysis 
values. 

The  initial  deviations  from  the  ideal  cylinder  are  ^11  =0.02  for 

the  upper,  0.05  for  the  middle,  and  0.12  for  the  lower  regions  of  the 
bandwidth.  These  represent  maximum  strength  to  classical  buckling 
strength  ratios  of  0.79,  0.6Q,  and  0.55  respectively.  For  comparison 
purposes,  these  imperfection  parameter  values  are  retained  through  the 
inelastic  analysis. 

INELASTIC  LOAD-SHORTENING  RELATION 


Addition  of  the  nonlinear  terms  associated  with  nonzero  values  of  the 
Ramberg-Osgood  parameter  K  to  the  linear  elastic  load-shortening 
relation  for  a  two-element  isotropic  circular  cylindrical  shell  re¬ 
sults  in  the  family  of  radius- to-thickness,  material-dependent  curves 
(solid  lines)  in  Figures  5  through  8.  Each  figure  represents  an 
imperfection-sensitive  series  of  curves  for  a  particular  material 
(aluminum  or  stainless  steel)  and  R/t  ratio  (200  or  100).  The 
predominant  failure  mode  remains  the  same,  although  the  maxima  8re 
significantly  affected  relative  to  the  aluminum  cylinder  of  R/t  =  200. 

The  major  effects  exemplified  in  the  two  choices  of  R/t  and  material 
are  the  transition  from  purely  elastic  to  elasto-plastic  behavior  with 
R/t  decrease  and  the  significant  reduction  in  maximum  strength  with 
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Figure  9.  Experimental  Buckling  Loads  for  Axially  Compressed 
Circular  Cylindrical  Shells. 


decreasing  values  of  N  in  the  uniaxial  constitutive  law.  In  Figure  5, 

2024-T3  aluminum  with  an  R/t  =  200  can  be  considered  to  act  within  99$ 

of  the  elastic  values  at  all  times.  The  R/t  dependency  in  Figure  6 

appears  as  a  direct  result  of  inelastic  effects.  The  13$  reduction 

in  maximum  strength  for  R/t  =  100  and  ^11  =0.02  as  opposed  to 

£  ® 

R/t  =  200  and  5llo  =  0.02  in  Figure  5  results  from  the  occurrence 
of  local  stresses  exceeding  the  proportional  limit  stress  of  the  ma¬ 
terial  and  the  degree  of  material  nonlinearity.  In  Figure  8  and  9> 
the  same  R/t  -ratio  effects  can  be  observed,  but  with  as  great  as  a 
44$  reduction  in  maximum  strength  when  the  material  is  stainless  steel 
and  ^11q  remains  equal  to  0.02  .  The  highly  nonlinear  nature  (N  =  3) 
of  the  uniaxial  constitutive  law  produces  considerable  softening  in  both 
the  peaks  and  the  slopes  of  the  load-shortening  curves.  For  these 
relatively  low  values  of  R/t,  the  peak  stresses  are  significantly 
above  that  of  the  proportional  limit  although  the  lateral  wall  deflec¬ 
tions  are  still  of  the  order  of  the  wall  thickness.  The  analysis,  then, 
is  still  within  the  kinematic  limitations  of  von  Karman-Donnell  theory 
discussed  by  Mayers  and  Rehf ield. ^  The  relatively  low  R/t  ratios 
are  not  extreme,  since  practical  cylinders  (either  stiffened  or  sand¬ 
wich)  often  have  "effective"  radius-to-thickness  ratios  in  the  range 
50  -  250  .  Further,  the  0.2$  offset  yield  stress  is  far  from  exceeded. 

Figure  10  is  a  summary  curve  developed  from  the  static  maximum  strengths 
of  Figures  5  through  8  by  plotting  them  versus  initial  imperfection 
parameter  ^11q  .  The  linear- elastic  curve  represents  the  theoretical 

upper  bound  to  axial  compression  buckling  of  perfect  and  imperfect 
shells.  The  remaining  solid  curves  show  the  reduction  in  maximum 
strength  due  to  the  combined  effects  of  initial  imperfections  and  ma¬ 
terial  nonlinearity  as  a  function  of  R/t  ratio.  Included,  also,  are 
imperfection  "knockdown"  design  stresses  adjusted  for  plastic  effects 
as  recommended  in  Reference  30.  The  range  of  imperfections  shown  in 
Figure  10  represents  the  distribution  of  initial  imperfections  that  are 
normally  used  to  explain  the  high,  medium,  and  low  maximum  stresses 
observed  in  past  experimental  investigations  (see  Figure  7).  At  best, 
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Figure  10.  Prediction  of  Maximum  Strength  Versus  Initial- 
Imperfection  Parameter  for  2024-T3  Aluminum  and 
Stainless  Steel  (l/2  Hard)  Cylinders  (R/t  =  100 
and  200)  and  Comparison  With  a  Standard  Procedure 


the  "knockdown"  stresses  of  Reference  30  are  still  about  10$  below  the 
present  results  in  the  case  of  2024-T3  aluminum  and  R/t  >  200  .  For 
stainless  steel,  with  an  R/t  =  100,  they  are  about  40$  below  the 
present- analysis  predictions. 

These  curves  show  that  the  greatest  effect  in  maximum  strength  reduc¬ 
tion  may  well  be  the  initial  imperfection  sensitivity  propounded  in 
the  literature,  but  that  in  the  range  of  practical  shell  "effective" 
radius -to- thickness  ratios,  material  nonlinearity  can  significantly 
affect  maximum  strength  performance  at  stress  levels  well  _>elow  those 
defined  by  the  0.2$  offset  yield-stress  criterion.  Further,  the  effect 
of  material  nonlinearity  can  contribute  to  the  scatter  of  test  results 
at  a  given  ratio  of  R/t,  as  discussed  in  Reference  31* 

ELASTIC  DYNAMIC  PERTURBATION 


The  effect  of  a  dynamic  perturbation  on  the  static  maximum  strength  of 
a  uniformly  shortened  shell  is  assessed  by  incorporating  the  mechanism 
of  time -dependent,  initial-imperfection  growth  in  the  static  analysis 
described  in  Appendix  II.  With  inertia  effects  neglected,  this  analog 
is  seen  from  Figures  5  and  6  to  yield  a  first  approximation  to  the 
actual  effect  of  nonimpulsive  dynamic  disturbance  of  an  imperfect 
shell  under  constant-rate  end  shortening.  The  growth  of  the  initial 
imperfection  amplitudes  in  sinusoidal  fashion  with  time  is  related  to 
the  end-shortening  rate  and  creates  a  quasi-dynamic  perturbation  of 
the  static  problem  as  described  in  Appendix  IV.  In  view  of  the  sen¬ 
sitivity  of  the  elastic  buckling  process  for  thin  shells  to  small 
deviations  from  the  perfect  shape,  the  analog  approach  obviously  offers 
a  much  more  tractable  solution  of  the  kinematically  nonlinear  problem 
and,  as  is  discussed  in  the  subsequent  section,  forms  the  basis  for 
extension  to  the  both  kinematically  and  constitutive!}  nonlinear  problem. 
Since  maximum  loads  are  reached  very  rapidly  with  small  increases  in  im¬ 
perfection  amplitude,  it  is  not  necessary  to  consider  more  than  the 
first-quarter  period  of  sinusoidal  oscillation  following  application 
of  the  disturbance. 
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With  the  maximum  amplitude  of  dynamic  imperfection  growth  fixed  at 
e  £ 

11  +  20  =  0.0125,  static  linear  elastic  maximum  strength  load- 

°  t 

shortening  curves  are  constructed  for  Ml  =  0.02,  0.05,  and  0.12 
as  shown  in  Figure  5»  An  aluminum  cylinder  of  R/t  =  200  behaves  in 
linear  e’astic  fashion;  thus,  Figure  5  is  valid  for  all  linear-elastic 
cylinders  of  R/t  £  200  within  the  framework  of  von  Karman-Donnell 
shell  theory.  Concomitantly,  static  load -shortening  curves  are  also 
constructed  for  ^11  =  0.03,  0.0 6,  and  0.13,  respectively.  On 

the  upper  curve  in  each  set,  the  sinusoidal  imperfection  growth  is 
seen  to  commence  at  an  eR/t  value  corresponding  to  the  maximum 
value  of  crR/Et  on  the  lower  curve  in  each  set.  The  intermediate 
curves  represent  the  effect  of  dynamic  perturbation  of  the  initial 
imperfection  pattern.  The  peaks  of  the  intermediate  curves  vary  as 
a  function  of  the  rate  of  application  of  unit  end  shortening  and  the 
rate  of  amplitude  increase  in  the  dynamic  perturbation  of  the  initial 
imperfection  pattern.  For  incremental  changes  in  eR/t  of  0.05, 

0.025,  and  0.010,  respectively,  the  circular  frequencies  of  the 
sinusoidal  oscillation  are  cn  =  0.188,  0.377,  and  0.942  rad/sec, 
respectively.  These  frequencies  are  well  out  of  the  range  of  natural 
frequencies  of  even  completely  unloaded,  simply  supported,  long  shells 
of  practical  proportions. 

From  the  curves  of  Figure  5,  it  is  noted  that  more  rapid  application  of 

the  disturbance  than  that  corresponding  to  a  change  in  eR/t  of  0.01 

in  about  1.7  seconds  would  lead  to  an  immediate  maximum  in  crR/Et. 

For  the  rates  shown,  however,  the  shell  either  maintains  its  load 

or  continues  to  load  until  a  maximum  is  reached.  The  maximum  load 

is,  in  general,  bounded  by  the  static  analysis  maximum-strength  values. 

For  relatively  small  imperfections,  the  practical  case  of  interest,  it 

is  seen  that  the  quasi-dynamic-analysis  maximum  loads  are  appreciably 

reduced  from  those  of  the  static  analysis  and,  further,  that  the  maxiita 

occur  well  prior  to  the  amplitude  of  a  disturbance  reaching  its  fixed 

maxima  of  ^11  =  0.01  (^11  +  ^20  =  0.0125)  •  'For  such  imperfection 
o  o  o 

growth  rates,  this  point  is  reached  when  an  intermediate  curve  joins  the 


lower  static  curve  in  each  set.  For  small  imperfections,  this  region 
is  not  shown  because  the  von  Karman- Donne  11  theory  hypotheses  have 
extended  their  validity  (see  Reference  10). 

The  results  shown  in  Figure  5  for  linear-elastic  behavior  would  appear 
to  indicate  that  imperfection-sensitive  shells  are  further  sensitive 
to  slowly  applied  disturbances  of  the  initial  imperfection  pattern 
occurring  in  the  vicinity  of  the  static  maximum  load.  This  phenomenon 
cannot  be  overlooked  in  evaluating  static  buckling  criteria  for  shells. 
Dynamic  disturbances  are  always  present  in  flight  environments  and,  to 
a  lesser  degree,  in  the  laboratory.  The  rationalization  of  large  scat¬ 
ter  in  compressed  shell  test  data  to  the  presence  of  initial  imperfec¬ 
tions  (see,  for  example,  References  1,  2,  and  3)  and  the  obviously 
statistical  nature  of  static  initial  imperfections  when  viewed  in  the 
light  of  the  results  presented  here  suggest  that  dynamic  disturbances 
of  relatively  low  frequency  occurring  during  the  loading  history  can¬ 
not  be  Ignored. 

INELASTIC  DYNAMIC  PERTURBATION 


The  effect  of  nonimpulsive  dynamic  disturbance  of  the  static  loading 
process  for  linear  elastic  shells  has  been  discussed  in  the  previous 
section.  However,  of  interest  in  the  present  analysis  is  the  effect 
of  dynamic  disturbance  coupled  with  the  nonlinear  elastic  (or  inelastic 
as  defined  herein)  behavior  ! n  maximum  strength  analysis  of  compressed 
shells.  With  the  use  of  the  same  unit  end  shortening  and  imperfection 
pattern  growth  rates  of  the  linear  elastic  analysis  (Figure  5);  the 
inelastic  analysis  results  are  given  in  Figure  6  for  2024-T3  aluminum 
with  R/t  =  100  and  in  Figures  7  and  8  for  stainless  steel  (j~  hard) 
with  R/t  -  200  and  100,  respectively.  The  uniaxial  stress-strain 
curves  used  in  the  analyses  for  2024-T3  aluminum  alloy  and  stainless 
steel  (^  hard)  are  shown  in  Figures  3  and  4,  respectively.  The 
figures  show  the  same  basic  trends  as  those  observed  in  th  linear- 
elastic  analysis  with  respect  to  the  effects  of  slowly  .pplied  dynamic 
disturbances.  As  shown  originally  in  References  14  and  15  by  Mayers 
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and  Wesenberg,  the  effect  of  material  nonlinearity  in  materials  of 
markedly  different  values  of  K  and  N  is  quite  obvious.  The  ma¬ 
terial  nonlinearity  lowers  the  maximum  stress  parameter  relative  to 
the  linear-elastic  prediction,  with  the  stainless  steel  cylinders 
being  much  more  sensitive  to  stress-strain  curve  nonlinearity  than 
those  made  of  aluminum.  The  general  softening  effect  of  material 
nonlinearity  tends  to  decrease  the  influence  of  dynamic  disturbance 
relative  to  the  linear-elastic  results.  This  is  understandable  be¬ 
cause  of  the  lower  <rR/Et  values  involved  when  the  maxima  oc:ur. 
Another  point  of  importance  is  the  influence  of  elevated-temperature 
environments;  in  general,  material  nonlinearity  effects  would  be 
noticeable  at  higher  radius-to-thlckness  ratios. 
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CONCLUDING  REMARKS 


The  study  of  the  effects  of  small  nonimpulsive  dynamic  perturbations 
of  the  initial  imperfection  pattern  of  axially  compressed,  circular 
cylindrical  shells  has  been  undertaken  with  the  use  of  static  maximum 
strength  theory  and  analysis  and  a  suitable  analog  to  account  for  ini¬ 
tial  imperfection  growth  during  the  loading  process.  The  dependency 
of  the  analysis  on  nonlinear  material  properties  is  obtained  with  the 
use  of  stress-strain  curves  for  2024-T3  aluminum  and  stainless  steel 
hard).  A  family  of  load- shortening  curves  is  developed  for  a 
representative  range  of  static  initial  imperfections  relative  to 
each  material  and  two  radius-to-thickness  ratios:  100  and  200.  A 
similar  set  of  load-shortening  curves  has  been  generated  to  show 
the  effects  on  maximum  strength  of  dynamic  perturbations  in  the 
static  initial  imperrection  pattern  occurring  during  the  loading 
process. 

From  the  present  results  for  both  elastic  and  inelastic  behavior  in 
the  bending  and  buckling  of  initially  imperfect,  axially  compressed, 
circular  cylindrical  shells,  consistent  with  the  materials  and  geome¬ 
tries  studied,  it  can  be  concluded  that  slowly  applied  dynamic  per¬ 
turbations  can  be  still  another  phenomenon  contributing  to  the  reduction 
in  classical  buckling  strength  and  the  scatter  in  experimental  data  at 
a  given  value  of  radius-to-thickness  ratio,  especially  in  the  range  of 
"effectively"  thick  shells  (stiffened  and  sandwich  construction).  With 
most  investigators  pursuing  the  perplexing  behavior  of  compressed  shells 
from  the  point  of  view  of  linear-elastic  analysis,  static  imperfection 
sensitivity,  and  statistical  techniques  for  identifying  deleterious 
imperfection  patterns,  it  is  important  to  call  attention  to  the  real- 
world  effects  of  material  nonlinearity  and  dynamic  disturbances.  The 
latter  phenomenon  bears  further  investigation  on  the  basis  of  a  true 
dynamic  inalysis  (inertia  effects  included)  in  view  of  the  results 
presented  herein  developed  from  an  analog. 
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APPENDIX  I 


EQUATIONS  OF  MOTION  AND  BOUNDARY  CONDITIONS  DERIVED 
_ FROM  MODIFIED-REISSNER  FUNCTIONAL _ 


The  Reissner  functional  for  prescribed  end  shortening,  as  contrasted 
with  end  loading,  is  simply 

ll"  =  III  (cr  e  +  a  e  +  t  7  -  f')  dV  (24) 
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The  von  Karman-Donnell  strain-displacement  relations,  modified  for  the 
two-element  cylinder  of  Figure  1  and  including  initial  radial  imperfec¬ 
tions,  are 
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The  stresses  in  the  top  and  bottom  faces  of  the  two-element  cylinder  in 
terms  of  the  midsurface  and  bending  components  are 
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After  integration  through  the  thickness  of  the  two-element  cross-section, 
the  Reissner  functional,  Equation  (2k),  becomes 
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Substitution  of  Equations  (25)  and  (26)  into  Equation  (27)  yields 
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The  average  stresses  and  bending  moments  are  defined  as 
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After  substitution  for  the  average  stresses  and  bending  moments.  Equa¬ 
tions  (29)  and  (30),  and  the  strains,  Equation  (25),  the  Reissner 
functional  becomes 
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The  kinetic  energy  of  a  freely  vibrating  two-element  circular  cylin¬ 
drical  shell  is  defined  as 
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after  integration  through  the  thickness  of  the  cross  section. 


With  a  Lagrangian  written  in  terms  of  the  Reissner  functional  l f  and 
the  kinetic  energy  f  as  (U "  -  T"  )  ,  the  enforcement  of  Hamilton's 
principle  as  denoted  by 


Km 

8  j  (U"  -  T")  dt 


(54) 


can  be  undertaken  with  simultaneous  variation  of  the  Equation  (34)  with 
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the  resulting  variation  yields 
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With  suitable  integration  by  parts,  Equation  (}6)  becomes 
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For  Equation  (57)  to  vanish  simultaneously  for  arbitrary  variations  in 
the  states  of  generalized  stress  and  displacement,  each  of  the  terms 
must  vanish  independently.  The  resulting  equations  of  motion  and  con¬ 
stitutive  relations  are  as  follows: 

Equations  of  motion: 
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Substitution  of  Equations  (38)  and  (39)  into  Equation  (4o)  simplifies 
the  lateral  equilibrium  equation  to  the  form 
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Stress-displacement  relations: 
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Moment-curvature  relations: 
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The  boundary  condition  combinations  involving  generalized  stress  and 
displacement  at  the  ends  of  the  shell,  x  =  0  and  x  =  L  ,  are 
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With  inertia  effects  neglected,  enforcement  of  Hooke's  law,  and 
2t^  -*  t  (the  thickness  of  a  homogeneous,  isotropic  shell),  the  set 
of  equations  and  boundary  conditions  given  by  Equations  (38),  (39  )> 
and  (hi)  through  (51)  reduces  to  the  von  Karman-Donnell  shell  theory. 
With  enforcement  of  Hooke’s  law,  2t^  -•  t,  and  the  inertia  terius 
due  to  motion  in  the  axial  and  circumferential  directions  neglected 
in  Equation  (hi)  in  accordance  with  the  conclusion  of  Reissner,^ 
the  set  of  equations  is  that  governing  the  nonlinear  free  vibrations 
of  shallow  shells  given  in,  for  example,  Reference  20. 
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APPENDIX  II 

METHOD  OF  SOLUTION  FOR  TWO-ELEMENT, 
LINEAR-ELASTIC  CIRCULAR  CYLINDER 


Substitution  of  the  assumed  radial  and  midsurface  displacement  func¬ 
tions  (Equations  (18)  and  (19)),  the  assumed  midsurface  and  bending 
stress  distributions  (Equations  (20)  and  (21)),  and  the  prescribed 
initial  imperfection  pattern  (Equation  (22))  into  the  modified 
Hamiltonian  represented  by  Equation  (l6),  followed  by  integration 
over  the  volume  of  the  cylinder  material,  yields 
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where  (U*)^  ^  is  the  contribution  of  nonlinear  material  behavior 

in  the  Reissner  formulation.  It  is  given  by 
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where  °"e££  Is  defined  in  terms  of  the  membrane  and  bending  stresses 
through  Equations  (10)  and  (26). 


For  the  elastic  case,  K  =  0  and,  therefore,  (U//).  ,  vanishes. 

In  the  resulting  set  of  equations,  the  15  generalized  stress  coefficients 
can  easily  be  expressed  explicitly  in  terms  of  the  free  average  stress 
and  displacement  parameters  through  enforcement  of  Hooke's  law.  The 
system,  as  a  result,  is  reducible  to  only  four  equations  in  terms  of 
the  four  independent  variables  ,  |2Q  ,  ,  and  aR/Et  and  may  be 

solved  numerically  for  prescribed  values  of  e  as  well  as  ^11q  ,  ^2C>o  , 
4  ,  and  t)  ,  the  initial  waveform  parameters;  a  relationship  between 
(oil/ Et)  and  (eR/t)  can  then  be  constructed  in  the  form  of  a 
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load-shortening  curve.  Although  ^20q  is  an  initial  waveform  parameter, 
in  view  of  the  justifications  reported  in  References  20,  26,  and  29, 
is  taken  equal  to  ^11q/4  throughout  the  present  analysis. 

With  the  inertia  effects  included,  all  arbitrary  coefficients  of  the 
assumed  displacement  and  stress  distributions  become  functions  of  time. 
Thus,  the  variational  process  in  such  a  case  leads  to  a  set  of  time- 
dependent  ordinary  differential  equations.  These  equations  together 
with  given  initial  and  loading  conditions  define  a  particular  problem. 
For  relatively  slow  end  shortening,  the  static  maximum  strength  solu¬ 
tion  is  assumed  to  be  little  affected  by  the  dynamics  of  the  buckling 
process  for  the  initially  imperfect  shell.  As  a  result,  this  type  of 
problem  is  not  of  interest  herein  as  is,  of  course,  the  impulsive 
loading  problem.  Of  concern  is  the  effect  of  a  dynamic  perturbation 
of  the  initial  imperfection  pattern  as  the  shell  is  shortened  into  the 
vicinity  of  the  end-shortening  magnitude  corresponding  to  maximum  load. 
The  effect  could  be  represented  by  retaining  inertia  effects  in  the 
problem  formulation,  adding  a  periodic  component  to  the  existing  ap¬ 
plied  constant-rate  unit  end  shortening  e  =  et  ,  and  specifying  the 
time  and  initial  conditions  relative  to  the  appearance  of  the  dis¬ 
turbance.  In  view  of  the  complexity  of  the  static  maximum- strength 
analysis  and  the  sensitivity  of  the  load-shortening  curve  to  initial 
imperfections  aid  nonlinear  material  behavior,  it  appears  reasonable 
to  consider  the  presence  of  a  small  dynamic  disturbance  by  neglecting 
inertia  effects  and  introducing  an  analog  to  estimate  the  first-order 
influence  of  such  a  disturbance  on  maximum  load.  The  analog  is  de¬ 
scribed  in  Appendix  IV. 

When  the  inertia  effects  are  neglected  and  the  analog  is  not  in  use, 
the  problem  solution  yields  identically  the  results  obtained  previously 
by  Mayers  ct  al.^’^  for  imperfect,  inelastic  (two-element  model)  and 
elastic  (conventional  model)  unstiffened  shells,  respectively.  When 
the  stress-strain  law  is  linear  and  imperfections  are  absent  (perfect 

shell),  the  results  of  References  14  and  28  are  identical  and  equivalent 

25 

to  those  first  obtained  by  Kempner. 
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REISSNER  FUNCTIONAL  FOR  THE  CIRCULAR  CYLINDRICAL 
SHELL  ~  THE  EQUIVALENT  STATIC  PROBLEM 


With  time  dependence  and  inelastic  effects  neglected  (K  =  0)  in 
Equation  (52)  of  Appendix  II,  the  Reissner  functional  reduces  to 
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This  functional  is  identical  to  that  developed  originally  in  Reference 

14  . 
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APPENDIX  IV 

DYNAMIC  ANALOG  CONSIDERATIONS 


The  use  of  an  analog  Co  assess  Che  effecCs  of  dynamic  perCurbaCions  on 
Che  maximum  sCrengCh  of  inicially  imperfect,  uniformly  shortened,  cir¬ 
cular  cylindrical  shells  is  demonsCraCed  in  Chis  appendix.  The  analog 
makes  use  of  Che  concepC  of  Cime-dependenC,  inicial-imperfeccion  growCh 
occurring  in  Che  buckling-sensitive  region  of  a  compressed  shell  Co 
represenC  Che  addicion  of  a  small,  nonimpulsive  dynamic  disCurbance 
Co  an  essenCially  sCaCic  loading  siCuaCion  for  a  flighC  vehicle  struc- 
Cural  shell.  The  effecCs  of  small  disCurbances  on  maximum  sCrengCh 
predicCions  for  specific  aluminum  and  sCainless  sCeel  cylinders  are 
shown  in  Figures  5  Chrough  8  and  summarized  in  Figure  10. 

As  noCed  from  EquaCion  (22),  Che  addicional,  Cime-dependenC  inicial 
imperfecCion  is 


This  perCurbaCion  in  Che  convenCional  (sCaCic)  maximum  sCrengCh  analysis 
provides  a  reasonable  firsC  approximaCion  to  the  occurrence  of  inertia 
effects  and  a  qualitative  appraisal  of  the  reduction  in  shell  maximum 
strength  to  be  expected  either  in  the  laboratory  or  flight  environment 
where,  at  least,  small  inertial  loads  can  easily  occur  in  a  statically 
loaded,  imperfect  shell  structure  loaded  in  the  vicinity  approaching 
the  static  stability  point. 

With  reference  to  an  end-shortening  controlled  laboratory  test  machine, 
the  unit-shortening  rate  e  is  assumed  to  be  3  x  10  ^  in. /in. /sec. 

Such  a  rate  would  cause  classical  buckling  of  a  perfect  shell  of 
R/t  =  200  to  occur  in  about  100  seconds;  the  rate  is  sufficiently 
8 low  as  to  preclude  any  considerations  of  buckling  under  relatively 
short-duration  loading. 
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For  initially  imperfect  shells,  defined  by  given  values  of  *11  (it 

t  e  O 

is  assumed  that  520q  =  S11QA  (see  Appendix  II)),  static  load¬ 
shortening  curves  are  developed  for  three  pairs  of  cases:  *11q  -  0.02 

and  0.03  ,  *11  =  0.05  and  0.06  ,  and  *11  =  0.12  and  0.13  •  As  in- 
’  o  ’  o  . 

dicated  in  "DISCUSSION  AND  RESULTS*’,  the  first  value  of  *11  in  each 
pairing  has  been  selected  to  correlate  maximum  strength  predictions 
with  the  bandwidth  of  experimental  data  at  R/ t  =  200  (see  Figure  9)- 
The  latter  value  of  *11  in  each  pairing  is  simply  a  static  refer¬ 
ence  bound  with  which  the  effects  of  varying  rate,  time-dependent  in¬ 
itial  imperfection  growths  can  be  compared.  The  onset  of  Initial 
imperfection  growth  is  assumed  to  occur  at  a  value  of  end  shortening 
parameter  eR/t  which  gives  the  stress  level  on  the  first  *11 
curve  corresponding  to  the  maximum  stress  level  on  the  second 
*11q  curve  in  a  set  of  pairings.  For  each  of  the  three  sets,  three 
circular  frequencies  of  oscillation  ire  selected  consistent  with  the 
unit  shortening  rate  of  3  X  10  ^  in. /in. /sec.  These  are  illustrated 
in  Figure  11(a),  where  the  end  shortening  parameter  eR/t  is  plotted 
versus  the  nondimens ional  time  parameter  u>t  .  The  growth  of  the 
total  initial  Imperfection  parameter  (*11q  +  *20q)  to  a  value  of 
0.0125  in  a  quarter-period  of  the  sinusoidal  oscillation  is  shown  in 
Figure  11(b),  also  as  a  function  of  cot.  Only  a  quarter-period  need  be 
considered,  since  the  maximum  load  of  the  dynamically  perturbed  shell 
occurs  within  this  range.  A  cross  plot  of  imperfection  growth  to 
maximum  amplitude  (*11q  +  *20q)  as  a  function  of  end- shortening 
parameter  eR/t  obtained  from  Figure  )+  is  presented  in  Figure  12. 
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